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Recently, there are two hints arising from physics beyond the standard model. One is a pos- 
sible energy loss mechanism due to emission of very weakly interacting light particles from white 
dwarf stars, with a coupling strength ~ 0.7 x 10~ 13 , and another is the high energy positrons ob- 
served by the PAMELA satellite experiment. We construct a supersymmetric flipped-SU(5) model, 
SU(5)xU(l)x with appropriate additional symmetries, [U(l)if] gaU ge x [U(l)_RxU(l)r] g iobai x Z2, 
such that these are explained by a very light electrophilic axion of mass 0.5 meV from the spon- 
^\ • taneously broken U(l)r and two component cold dark matters from Z2 parity. We show that in 

f*"^ ' the flipped-SU(5) there exists a basic mechanism for allowing excess positrons through the charged 

SU(5) singlet leptons, but not allowing anti-proton excess due to the absence of the SU(5) singlet 
C^l ■ quarks. We show the discovery potential of the charged SU(5) singlet E at the LHC experiments 

<— 1 ' by observing the electron and positron spectrum. With these symmetries, we also comment on the 

I mass hierarchy between the top and bottom quarks. 
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CLf I. INTRODUCTION 

^ n 

D t The most looked-for particles on Earth at the present time are axions [1[ and weakly interacting massive parti- 
cles(WIMPs) 0. Recently, hints on these particles have been reported from outer space sources. The 10 year old 
' '' INTEGRAL data gives the 511 kcV line. The white dwarf cooling has been suggested by the emission of very 
weakly interacting light particles (VWLP) (mass less than ~ 1 cV) with a very small coupling strength to electron 
£> ' (~ 0.7 x 10~ 13 ) [J]. At present, this white dwarf bound can be considered just as an upper bound, but in this 
1 \ paper we adopt a bold assumption that their best fit corresponds to the existence of a VWLP. More recently, the 
remarkable observation of high energy positron excess in the satellite PAMELA data [| hints heavy (100 GeV - 10 
TeV) cold dark matter(CDM) particles 0,0,11. Since the 511 keV line can be explained by the astrophysical origin, 
or by the particles of mass 0(1 — 10) MeV or by an excited state almost degenerate to the WIMP in this 
paper we attempt to understand the latter two of these observations from a grand unification (GUT) viewpoint of 
the flipped-SU(5) [II 03. 

The PAMELA group [5[ reports the high energy positron excess, above 10 GeV up to 60 GeV, at the level e + /(e + + 
e~) O(O.l), but has not reported any noticeable excess of anti-protons The charactersistic of the data between 
10-60 GeV is a slightly rising positron flux, which is inconsistent to the lightest supersymmetric particle(LSP) in the 
minimal supersymmetric standard model(MSSM) with supersymmetric(SUSY) particles at O(100 GeV) Thus, 
the MSSM with O(100 GeV) lightest neutralino DM component x( LN x) has to be modified if SUSY has to be kept 
for the gauge hierarchy solution. The minimal extension needs just three two component fermions Nr, Eh and E R 
Q, with the coupling 

W = e R E c R N R . (1) 

If one tries to introduce more fields around TeV scale, certainly the rising positron flux can be explained, but it is 
very difficult [H to explain nonobservation of antiproton flux. This suggests a different treatment of electron from 



quarks. Indeed, the flipped-SU(5) GUT SU(5)xU(l)x treats charged leptons differently from quarks [llj such that 
the charged lepton is an SU(5)xU(l)x singlet, 
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The specific form of the coupling |T|) is related to the problems of the Yukawa couplings. In this regard, we note the 
five old but fundamental problems on the masses of the standard model fermions. Firstly, there is a hierarchy of the 
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top quark mass m t ~ fv/\/2 ~ 170 GeV at the electroweak scale of v ~ 250 GeV, and the rest of the standard model 
particle masses are much smaller than the top quark mass, m < nib — 4.5 GeV <C nit- Second, there is a sizable 
Cabibbo mixing, especially between the first and the second family quarks. Third, the quark mass ratio is reversed 
in the first family compared to those of the second and the third family quarks, i.e. m u /md < 1 while m c /m s > 1 
and mt/rrib > 1. Even though it is not proper to define an inverted mass pattern for leptons, it looks like that the 
electron mass is also an inverted pattern because m e /m^ ~ m„/4m c ~ 0(1/200). Fourth, there is a big hierarchy 
of singlet neutrino mass of m Uo ~ O(10 14 ) GeV for the seesaw mechanism [Tj] and the electroweak scale of v ~ 250 
GeV. Finally, we point out that the SM neutrino mass pattern shows a tri-bi maximal mixing pattern [l5| . Except 
the tri-bi maximal mixing, in this paper we try to understand the remaining four problems in a flippcd-SU(5) GUT. 

To have a very light axion implied by the white dwarf data, we introduce a Peccei-Quinn (PQ) symmetry U(l)r 
so that the symmetry we consider is extended to SU(5)xU(l)x times U(l)r- To have the inverted mass pattern of 
the quark mass ratio in the first family, we assign different T charges for the first family members and the second and 
the third family members. But, to have a large Cabibbo mixing, we require Q em = § quarks or Q cm = — | quarks 
mix fully, restricted only by the strength of the corresponding Yukawa couplings but without any restriction from the 
symmetry arguments. In the flipped-SU(5), since both the quark singlet cLr and the quark doublet q R appear in 10i, 
viz. Eq. it is easy to mix Q em = — ^ quarks fully without restrictions from the symmetry arguments. So, we 

assign the same symmetry charges for three members 10^ (i — 1, 2, 3). Then, Q em = — | quarks can be mixed fully as 

desired. But the charges of 5^ 3 and l^ can be different for different i to allow the inverted mass pattern of the quark 
mass ratio for the first family. For the singlet neutrino mass problem, we need a large VEV, above the axion window 
of 10 9 — 10 12 GeV. If there is no other gauge symmetry, the PQ symmetry may be broken by this large VEV. Thus, 
we introduce another U(l) gauge symmetry which will be called a horizontal gauge symmetry U(l)#, distinguishing 
the first family and the second and third families. Finally, to keep Er — E R pair at the electroweak scale and to have 
proton stability and two dark matter components Q, we introduce the R symmetry and Z2 symmetry (the i?-parity). 
Thus, the symmetry we consider is 



In this model, the renormalizable couplings are gauge couplings and the top quark Yukawa coupling. All the other 
Yukawa couplings are non-renormalizablc. The model we introduce docs not have the C3 term of Ref. [l|, where the C3 
term is the axion and QCD anomaly coupling, and the PQ charges of Higgs doublets vanish. Therefore, the very light 
axion we obtain is not a proper Kim-Shifman-Vainstein-Zakharov (KSVZ) axion [l6|, or the Dine-Fischler-Srednicki- 
Zhitnitskii (DFSZ) axion [17| . but can be called a 'variant very light axion'. The electron-axion coupling is known 
to be small for the KSVZ axion, as calculated in Ref. [l8j|. In our case, the axion arises from the ci term of Ref. 
[l| , where the C2 term is the light quark mass term, through a non-renormalizable interaction, and the electron-axion 
coupling turns out to be large. Since this 'very light variant axion' is required to have a large electron coupling, it 
will be called a very light electrophilic axion. 

In Sec. UH we briefly review the white dwarf bound In Sec. IIII1 we present a flippcd-SU(5) model, satisfying 
the criteria presented above. In Sec. IIV1 we discuss the SUSY phenomenology for the model presented in Sec. IIIII 
Sec. |V] is a conclusion. 



White dwarfs can give us hints about their last stage evolution. White dwarfs, with the mass range less than 
I.O8M0, are the fossil stars after consuming all their nuclear fuel. Here the electron degeneracy pressure is the one 
supporting the gravitational collapse, and the surface temperature ranges down from ten times the solar surface 
temperature. The temperature difference between the core and the surface is not greater than 100 K, and hence if 
the VWLP is produced at the core then its mass cannot be much greater than 1 eV. Since the white dwarf has a 
relatively simple structure and its evolution mechanism is just based on the cooling, its luminosity function has been 
rather accurately predicted. Thus, if there is a novel energy loss mechanism by a VWLP like axion, it can affect the 
resulting luminosity function of the white dwarf which can be used to constrain the properties of that particle or give 
a clue for its existence. 

The most strong bound on VWLP comes from the existence of the sharp cutoff at the faint part in the luminosity 
function (n) versus the bolometric magnitude (Mboi) plot. The sharp cutoff is due to the finite age of our Galaxy. 
The luminosity function is given by Q 



SU{5) x U{l)x x U(1) H x U(1) R x U(l) r x Z 2 . 



(3) 



II. WHITE DWARF BOUND 
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where 4>{M) is the initial mass function, ip is the birth rate of the progenitor star, Mboi is the bolometric magnitude, 
M is the mass of the progenitor star and t coo i = dt/dM^oi is the characteristic time scale of the white dwarf cooling. 
Although the luminosity function depends on the birth rate which is not yet exactly known, the bright part of the 
curve is insensitive to it up to the overall normalization. Thus, one can exploit the bright part to constrain the 
hypothetical cooling mechanism due to the VWLP. 

In this regard, recently the possible energy loss rate has been computed [If. If a very light boson with mass less 
than 1 eV is produced inside white dwarfs, its coupling to electron is estimated as 0.7 x 10 -13 . Any boson with this 
coupling strength to electron can explain the white dwarf evolution a la Ref. [J , using Ref. (l9| . 

Since axions are well-motivated light pseudoscalar bosons, Ref. 0] tries to interpret the coupling 0.7x 10~ 13 in terms 
of the DFSZ model [l7|- In the interior of the white dwarf, the axion bremstralung coming from the axion-electron 
Yukawa coupling, which is significant in the DFSZ model but absent in the hadronic axion model, is dominant over 
the Primakoff process. The energy loss rate per mass with the crystalization effect is given by 

-l.OSxlO-^^)^) 4 ^^) ergg-s- (5) 

where F(p, T) is the correction factor due to the crystalization effect and g aee is the Yukawa coupling between the 
electron and the axion. In the DFSZ model, it corresponds to 2.8 x 10~ n sin 2 j3(m a /leV) with the domain wall 
number 6. From this study, Isern et al. [1] point out that the small deviations of the white dwarf luminosity function 
from the standard cooling model can constrain the axion mass m a sin 2 (3 < 10 meV, and even presented the best fit of 
the luminocity function with the axion emission for m a sin 2 f3 ~ 5 meV. In the DFSZ model, this white dwarf bound 
on the axion mass corresponds to the axion decay constant F a > 0.6 x 10 9 sin 2 f3 GeV and the best fit corresponds to 
F ~ 1.2 x 10 9 sin 2 (3 GeV. Below, we try to construct a flippcd-SU(5) to accomodate the best fit of Ref. 0]. 

The electron axion coupling can be written as (m e T(e) / F)eij^ea, where F = N^wFa and T(e) is the PQ charge of 
e. Since the axion mass is given by 0.6 eV(10 7 GeV)/F a , the best fit of Isern et al. [4j in the DFSZ model interpreted 
as the electron-axion coupling in any axion model, 



Best fit 



m e r(e) 



F a 



m P 



0.72 x 10 10 GeV 



0.7 x 10 13 , in any axion model. 



(0) 



If the domain wall number is large and the PQ charge of electron is 1, then we need to have a smaller F a than 10 10 
GeV since NowF a is required to be 0.72 x 10 10 GeV. Generally, the domain wall number is counted as h for any 
quark carrying the unit PQ charge. Therefore, if only one Dirac quark carries the PQ charge 1 so that each chirality 
contributes one unit of the PQ charge, then the domain wall number is 1. In this way [J], the DFSZ model gives 
Ndw = 6 since there are six quarks, and the electron-axion coupling is m e /(6 x 1.2 x 10 9 GeV) with F a = 1.2 x 10 9 
GeV. 

There is an interesting case that only one Weyl quark carries the PQ charge 1. In this case, our formula for the 
axion-electron coupling must use Ndw — j- 1 In such a model Eq. ([6]) implies F a ~ 1.4 X 10 10 GeV. Below we 
concentrate on the latter possibility in a flipped SU(5) model, obtaining Eq. ^ for the very light electrophilic axion. 
This moves F a to the supergravity favored scale [20j even with the large electron-axion coupling envisioned in Ref. 



III. A FLIPPED SU(5) MODEL 

SO(10) GUT groups cannot realize our idea of Ref. 0. The SO(10) GUT with the SU(5) GUT assignment 
[2ll ] docs not introduce our SU(5) singlet E and E c fields since the electromagnetic charge operator Q em is fully 
embedded in SU(5)gg- The flipped-SU(5) assignment of hypercharge in the SO(10) subgroup [ll|, namely 16 = 
10i + 5_3 + I5, allows the SU(5) singlet (in I5) carrying the electromagnetic charge, and en does not belong to 
10i or 5_3. Supersymmetric (SUSY) generalization of the flippcd-SU(5) chain from SO (10), including the vacuum 
expectation value of Higgs field not belonging to the origin of the weight diagram to reduce the rank of the gauge group, 



has been discussed in 12j. A simple choice of such symmetry breaking with 10h and 10# from string construction 



was shown to be possible in the fermionic construction [22j and in the orbifold compactification [23j. Also, in the 



1 There cannot be a half integer domain wall number. It simply means that with a proper definition, the electron PQ charge is 2 compared 
to the quark charge 1. See Eq. H71 . 
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SUSY flippcd-SU(5), the doublet-triplet splitting problem is solved elegantly as shown in the fermionic construction 
[24| and in the orbifold compactification [23|. 

Here, we do not consider the SO(10) breaking chain to an SU(5)xU(l)x gauge group. The SU(5)xU(l)x gauge 
group is considered as an ultraviolet completion gauge group, possibly arising from a compactification of string theory. 
As in Ref . @] , we use the right-handed chiral fields so that the electron singlet is represented as en instead of e c L and 
the quark doublet is represented as q c R instead of qh- This right-handed chiral field notation is simpler to discuss the 
interaction e R E R N R which we try to introduce for the PAMELA data. 

We use the following right-handed notation for representations of the flippcd-SU(5). The electroweak hypercharge 

is 

r-\r,-\x (7) 

where the SM singlet generator embedded in SU(5) is 

The representations @ are consistent with this definition. 

The matter field 10i contains the quarks doublet (3*, 2)_!/ 6 , the quark singlet (3, l)-i/3 and the neutral lepton 
singlet (l,l)o. 5_3 contains the quark singlet (3, 1)2/3 and the anti-lepton doublet (1, 2) 1 , 2 - (1,1)5 is the charged 
lepton of the unit charge 1. Therefore, it is useful to remember that singlets da, ur and en are members of 10i, 
5_3 and (1, 1)5, respectively 



A. Breaking nipped SU(5) 

The singlet neutral direction (i/q and vq directions) of 10i and 10 1 can be assigned to GUT scale vacuum expecta- 
tion values(VEVs), breaking SU(5) xU(l)x down to SU(3)x SU(2)xU(l)y. For the spontaneous symmetry breaking 
of the flipped SU(5), thus, we need a vectorlikc pair, IOh and 10//, which contain the SM singlet components, 
and vq. The following superpotential, 

W H = ^(10 ff T0 H - Mq) (9) 

allow 10// and 10// to develop VEVs along the direction of Vq and Vq. Here the superfield Z is a singlet of the flipped 
SU(5). Since the D flat direction is along (10//) = (10//), thus, at the SUSY minimum we have 

(10 H ) = (T0 H ) = M G . (10) 

The phase modes of q H (SU(2) doublet) and qn (SU(2) doublet) in 10// and 10// correspond to the Nambu-Goldstone 
bosons and the real modes of q H and qn become heavy as the Higgs boson becomes heavy in the Higgs mechanism in 
the standard model. The Nambu-Goldstone bosons are absorbed by the heavy gauge sectors when the flipped-SU(5) 
is broken. However, dn (SU(2) singlet) and dn (SU(2) singlet) in 10// and 10// need to be made somehow heavy. 
The quantum numbers of Z. 10//, 10//, and some other supcrficlds needed later arc displayed in Tablc[I] 
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TABLE I: Vectorlike fields including the Higgs fields. Z and $ are the flipped-SU(5) singlet fields. 

The representation 5_2 contains Q em = ^ quark D and Higgs doublet Hd with Y = i. The representation 52 
contains Q om = — -| quark D and Higgs doublet H u with Y = — |. In order to decouple the modes D, D in 5_2, 5 2 , 
and also d//, dn in 10//, 10// from the low energy field spectrum, we consider the following superpotential, 

W D/T = -jrjr- (10j/10//5 h + T0ifl0K5 fc ). (11) 
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We assume that the singlet <£>, which carries the charge of the gauged U(1)jj, get a VEV of order Mq. We will provide 
later an explanation on how <E> can get the VEV. Because of the VEVs (IQh) and (10 h) in the vp and Vq directions, 
D in 5_2 and du in 10h, and also D in 52 and dn in 10# obtain superheavy Dirac masses p3l. |24||. 

B. Yukawa couplings and family hierarchy 

The top quark mass is close to the SU(2)xU(l)y breaking scale and hence t is required to obtain mass by the 
renormalizablc terms. Since the b and r masses are O(10 -2 ) of the top quark mass, it is reasonable to require for 
them to appear by nonrenormalizable terms. Thus, we allow the renormalizablc supcrpotcntial term for the top quark 
and nonrenormalizable terms for the bottom quark and the tau lepton masses due to the gauge, global and discrete 
symmetries we introduce. For the second family members, the third family pattern is repeated but we allow 0(1O -2 ) 
couplings compared to the third family couplings. However, we treat the first family quite differently from the second 
and third families because the mass pattern is inverted for the up quark, i.e. m u /rrid < 1 while m c /m s > 1 and 
mt/rrib > 1. We interpret the electron mass also by the inverted pattern because m e jm il ~ m u /4m c ~ 0(1/200). We 
investigate this kind of mass patterns by introducing the PQ symmetry U(l)r and the horizontal gauge symmetry 
U(l) ff . 

As explained above, the 5_3 representation of Eq. in the flipped SU(5) accommodates a quark singlet ur 
(Qem = §) [instead of dn as in SU(5)gg] and a lepton doublet Ir (Qem = \)- As a result, considering only the 
SU(5)xU(l)x symmetry, u type quarks and Dirac neutrino Yukawa couplings can be obtained from 10i5_35^ with 
the dictated relation = muv in the flipped-SU(5) [rather than mj = m e as in the SU(5)gg]- Since the singlet 
lepton ep remains also as a singlet of the flippcd-SU(5) model unlike in the SU(5)gg> the Yukawa couplings for the 
charged leptons is given by 1^5— ^5h- The d type quarks' Yukawa couplings come from the couplings of the type 
lOilOiS^ in the fhpped-SU(5). 
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TABLE II: The matter fields and S. S is an SU(5) singlet field whose VEV is breaking the PQ symmetry spontaneously. 

Therefore, we assign the Abelian charges given in Tables U and [TTJ The Z2 symmetry distinguishes the matter lO' 2,3 - 1 
from the GUT breaking Higgs 10 H - With the U(l)j? charges of Tables U and M there is no U(1)#-SU(5)-SU(5) 
anomaly. Introducing exotic singlets with Q cm = — -|, £ and C in Table [I] the mixed anomalies of U(l)jy-U(l)|- 
and U(l) 2 ? -U(l)x are made to vanish. A possible TrX anomaly is made to vanish by introducing the H = and 
X = -1 singlets in Table [ml 

The Yukawa couplings for the second and third families of the MSSM matter fields, consistent with the symmetries 
of Tables U and HH are 

^(2,3) = JL 10 (2,3) 10 (2,3) Bh + io^s^s, + -^-1^5 {2 - 3) 5 h . (12) 
Mp Mp 

For the third family, therefore we qualitatively explain why the b quark mass is much smaller than the t quark mass, 
depending on the factor ($) /Mp. The mystery of the quark masses is not in the mystery of the heaviness of top quark, 
but resides in the smallncss of the masses of the first family members compared to the electroweak scale. Therefore, 
for the first family, we assign quantum numbers differently from those of the third family members. One simple choice 
is that the d quark mass is still made small by the choice of Yukawa couplings but the up quark and electron masses 
arise from different type of couplings. This choice seems reasonable since the up quark mass is somehow inverted from 
the second and the third generation pattern, m 2 /3 > m_i/ 3 . The electron mass is much smaller than the up quark 
mass and hence can be thought as the inverted pattern. Therefore, let us choose the quantum numbers of the first 
family members as shown in Table HT1 

In Table HH we also show the singlet S whose VEV breaks the PQ symmetry. The heavy fields 10', 10 ,1' and 1 of 
Table U will be integrated over to obtain the needed non-renormalizablc interactions. Note that for the d type quarks 
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the three flavors {dp, sp, bp} and {d R , s c R , b R } can be fully mixed among themselves, since all the quantum numbers 
of lCr > arc the same as those of 10 (2 ' 3) . This full mixing between Qom = — 1/3 quarks enables one to generate a 
large Cabibbo angle, even though Q C m = 2/3 quarks are not mixed fully. 

We assume that the VEV of the singlet S is around the intermediate scale 10 10 GeV. We will explain later how S de- 
velops such a VEV. The allowed nonrcnormalizable Yukawa couplings ((S) /AIp)W^5 {1) 5 h and ((5$) /Mj,) lW5 (1) 5 h 
by the gravity mediation gives, therefore, too small u quark and electron masses. Thus, we introduce vectorlikc pairs 
of {10', 10 } and {l',l } with proper quantum numbers given in Table [l] and their heavy mass terms. Then, the 
Yukawa couplings relevant to the u quark mass are 

W {u) = <M0 (1) T0~' + M'lO'TU' + 10'5 (1) 5 ft , (13) 

where the mass parameter M' can be smaller than Alp. After integrating out the heavy {10', 10 } pair, one can 
obtain the u quark Yukawa coupling with a desired size, 

W (u) = M 10 (i) 5(1)5^ (14) 
M' 

For the first family up-type quark, we have an extra factor ($) jAV compared to the second and third family up- 
type quarks. So, with an O(10 -1 ) smaller Yukawa coupling compared to the charm quark Yukawa coupling and 
(<S>)/M' ~ 0(1O -2 ), m u can be made smaller than m c by a factor 0(1O -3 ). 
Similarly, we have the following interactions for electron, 

W (e) = S1 (l)f + M'^rj' + * 1 /g< 1 ) 5fc< (l 5 ) 

Alp 

After integrating out 1' and 1 , the effective electron Yukawa coupling becomes 

(e) = _JS$)_ (1) ( } 
M"M P h v ' 

with (S$)/M"M P ~ 1CT 6 . 

This model does not have the C3 term of Rcf. [l[ and the PQ charges of Higgs doublets vanish. Therefore, the very 
light axion we obtain is not a proper KSVZ or DFSZ axion. Out of 10', lo', and 5^ , carrying PQ charges, disregard 
the vectorlike 10' and lo'. Then, there remains a half quark in 5 carrying the PQ charge since the other half in 
10 (1) does not. carry the PQ charge. After integrating out. heavy fields, this model therefore has the c 2 term of Ref. 
[l[ and is a variant very light axion as some electroweak scale axions were called variant axions (25| . The mass term 
related to C2 is a non-renormalizable interaction, and the electron-axion coupling will be shown to be large. The 
resulting c 2 type coupling is for the up quark mass term, m u ui,Upe %a ^ F + h.c. pj. Expanding S in terms of the axion 
field, S = (F + p)e 2la / F with T(S) = 2, Eq. (M]) gives the following axion-electron coupling, 



(\S\*H d ) 



M"M P 

= — m e ee + ^ m e^;^ evy^e a (17) 

In this model F a — F since effectively only one Weyl quark in 5 carries one unit of the PQ charge. This model 
gives effectively Np,w = 1/2 since only a half quark carries the PQ charge. See footnote [T] Note that a more general 
electrophilic axion can result with the PQ charge carrying heavy quarks giving C3 = iVh.q. and c 2 = — C3 ± i. 



C. Masses of singlet neutrinos 

For the seesaw generation of neutrino masses, it is necessary to have (i) the Dirac masses between the SM neutrinos 
and the singlet neutrinos at the electroweak scale and (ii) the heavy Majorana masses of C(10 14 ) GeV for singlet 
neutrinos. In string compactification to the MSSM, these two conditions are achievable in Z 3 ,Z 12 _/ and Z e _jj 
orbifold compactificaton [261 ] . 
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In the presence of the U(l)ij symmetry in the flipped-SU(5), it is not straightforward to write down the Majorana 
neutrino mass term of O(10 14 GeV) for the singlet neutrino because it is embedded in 10x- So, the R symmetry must 
be broken nontrivially by the VEVs of i?-charge carrying scalar fields [13]. We consider the following interactions: 

W (u) = — 10 (i) T0 H ** + -r^$$'* 2 , (18) 
Alp Mp 

where the quantum numbers of <I>, and \I/ are shown in Table Hill We introduced 3>" in Table [TTT1 to make TiH = 0. 
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TABLE III: The Q cm = -| singlet 77, and neutral singlets carrying nonzero U(l)ir or PQ charges. We need four 77s and one $" 
fields to make TrX = and TrH — 0. Fields in the brackets are already listed in Tables HI and HTl 

We assume that both <I> and <£>' get VEVs of order Mgut- Then, W achieves a heavy mass. After integrating out 
"J, one can get the effective Majorana neutrino mass term, 

2 

M = g)gOg) 10 (0 10 0-), (19) 
($')Mp 

Hence, if (<f>) ~ ($') and the dimensionless Yukawa coupling is of order unity, we get a Majorana neutrino mass of 
O(10^ GeV) as desired. 

£, £ and two of 77s can achieve masses by ($). To give masses to {£, C, ^(s)) ? 7(4)}j we need more fields which we omit 
here. We note that these neutral singlets can be almost massless since their decoupling temperature is expected to 
be very large. 



D. The PQ symmetry breaking scale 

Now let us discuss how $ and S can achieve the VEVs. Tables HI HTl and ITTTI list the fields carrying nonzero U(l)# 
charges. The relevant D term potential of U(l)# is 

v H = ^f\m 2 + \^'\ 2 -i*i 2 -e| 2 , (20) 

where £ (~ Mq) denotes the Fayet-Iliopoulos term [28[, and we have neglected the small contributions, |5/j| 2 , 15/J 2 , 
etc. Considering the soft mass terms, the minimum of the potential would be on |$| 2 + |<f>'| 2 = £. In addition, 
consideration on the higher order term in the Kahlcr potential, K D jp- |<E>| |$'| , and its supergravity effect on the 

potential would determine the minimum at ($) = ($') = (In supergravity, the Fayet-Iliopoulos term needs to 

be replaced by a proper modulus [29|.) 

Now, S can develop the intermediate scale VEV (~ 10 10 GeV) by the supergravity couplings and the softly broken 
SUSY effects, 

W PQ = j^S 2 S 2 . (21) 



The soft SUSY breaking mass terms, the A-term corresponding to Eq. (|2Tl) . and the F term contribution by Eq. ^1 
can make it possible to develop their VEVs at the minimum of the scalar potential, 



(S) ~ (S) ~ yJm 3/2 M P 



m-3/2 M P 



1/2 



10l ° GeV (ll^-H^) • W 
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At this scale, the PQ symmetry is completely broken. The gravitino mass of O(10 2 GeV) leads to the axion scale of 
O(10 10 GeV) we presented in Sec. M 



E. Two dark matter components 

In Ref. Q, two DM components are introduced. In fact, there arc two possible ultraviolet completions (items a 
and b below) with N R of Eq. (JTJ) . These are two high energy completion from the LNx and N R of of our model. 
These are based on how N R obtains mass. 

a. The needed matter fields in this case are SU(3) x SU{2) singlets E R , E r , N r , and N R . In the flipped-SU(5), 
they can appear as SU(5) singlets because the weak hypercharge depends on the SU(5) singlet generator X also, 
as shown in Eq. ([7|). The needed couplings are 0] 

W DM = e R E c R N R + N%, (23) 
where e R = I5 , and hence we assign the charges shown in Table HVl The R symmetry forbids superheavy 



Fields 


Er 


Er 


N R 


N R 


X 


-5 


5 








H 














R 


4 
3 


8 
3 


2 
3 


10 
3 


z 2 






+ 


+ 


r 


1 


-1 









TABLE IV: The SU(5) singlets for the coupling Q. They do not carry U(1)h charges. 

mass terms of E R E R and N R as well as H u Hd- Their mass term would be possible through the Giudice-Masiero 
mechanism (30l |: 

W D m 3/2 E R E c R + m' 3/2 N R N c R + m% /2 H u H d , (24) 

where m 3 / 2 , Wgy 2 , and m'^, 2 are of order the gravitino mass scale. The masses of E and N are at the 0(100 GeV) 
scale. 

b. There is another method to give mass to N R . It is by the R charge assignment of N R be | so that the A^^ 
term of (|23p is not allowed. The relevant quantum numbers are given in Table IVl In this model, the following 



Fields 


Er 


Er 


Nr 


S' 


X 


-5 


5 








H 














R 





4 


2 


4 


Z 2 


+ 


+ 




+ 


r 


1 


-1 









TABLE V: The SU(5) singlets for the coupling Q. S' is superheavy, 
supergravity effect generates the following mass terms, 

J d A 9^- (\E R E R + X'N R + X"H u H d ) . (25) 

At low energy, there are three more two component fields, E Rl E Rl and N R . 

The R charges of a are different from those of Ref. because in the present flipped-SU(5) model the R charge of 
H u Hd is necessarily nonvanishing while in Q it is zero. In the present case, the number of fields carrying the R 
symmetry breaking VEV can be one: S'. 

In the phenomenology discussion of Sec. IIV1 the details of the ultraviolet completion does not matter because the 
main prodution we will consider is through the electroweak process pp — > E~ E + and the E decay via the interaction 
(UJ) which is common to Models a and b. 



9 



IV. PHENOMENOLOGY 
A. PAMELA high energy positrons without antiproton excess 

The model presented has the coupling (p} and explains the high energy positron excess of the PAMELA experiment 
Q, but does not introduce a urU r N r type and/or djiD R Nji type couplings with heavy quark U and D. If such 
couplings were allowed, two DM annihilations would give high energy antiproton excess. But these terms cannot be 
present. Firstly, these heavy quarks carrying the color charge cannot be flippcd-SU(5) singlet representations. They 
must appear in 10i or 5_3 or some other SU(5) representation with a proper X quantum number. Because these are 
heavy quarks, say Q where Q is U or D, they must be vectorlikc below the GUT scale, Qr + Q C R . Both Q and Q c 
cannot appear in 10 or 10 because 10 ■ 10 • N and 10 • 10 • N are not allowed by the SU(5) symmetry. Similarly, both 
Q and Q c cannot appear in 5 or 5 because 5 ■ 5 ■ N and 5 • 5 • N are not allowed by the SU(5) symmetry. So, if one 
of Q and Q c appears in 10 (or 5), the other must appear in 10 (or 5) to write down QQN coupling. But we need 
a coupling with a light quark q where q is u or d. So, at least we must keep ten more two component fields at the 
electroweak scale, e.g. 5 and 5. So, assuming the minimal extension with N, it is just E + E c . Allowing 5 and 5 in 
addition will produce antiproton excess too but then we introduce ten more chiral fields. Even though we introduce 
5' and 5' and a coupling qRQ c R N R , deisregarding the minimality, if toq is sufficiently heavy (toq > M x + tun — m q ) 
then antiproton excess is not expected. So, the flippcd-SU(5) has an intrinsic mechanism to allow the positron excess 
without antiproton excess by an appropriate mass parameters. 

Beyond this naive expectation, in fact we can show that the diagonalization of mass matrix almost forbids the 
antiproton excess in our flipped-SU(5) model. For the heavy vcctorlike pair, we assign the quantum numbers of 
right-handed top for 5' and their opposite quantum numbers for 5', 

5"'_ 3 : X = -3, H = 1, R = 0, Z 2 = -, T = 
5+3 : X = 3, H = -1, R = 0, Z 2 = -, T = . 

The mass matrix of Q cm = 2/3 quarks can be written as the following matrix. Here m and m' are at the electroweak 





5« 


5-3 


loi 1 ' 


em 


m! 


5 3 








scale, to is a free parameter, presumably around 100 GeV, and e is of order 10 5 as discussed in Eq. (fT4|) . Then, the 
physical right-handed up quark and the heavy quark are given by 

ur = cos 9 + sin# 5'_ 3 (26) 
U R = - sin 6 1$ + cos 9 ]>_ 3 (27) 

where tan 9 — em/m' and 5'_ 3 represent the Q C m = 2/3 quark in 5'_ 3 . Then, the coupling urU r Nr is suppressed by 
e, and the antiproton excess is not expected in the flipped-SU(5) model. 

B. Production of E in LHC experiments 

The model presented in Ref. [7[ has three two component fermions. Among these Q cm = — 1 lepton E can be 
produced at the LHC machine. The most distinguishable feature of the model presented here is leptonic (especially 
'electronic') property of the dark matter component. E + and E~ can be produced by proton-proton collision and 
eventually they decay to an electron-positron pair and missing energy by the interaction (ft}. Then, we compare it 
with the expected high energy e + e~ spectrum from the LHC machine. This signal distinguishes the model of Ref. 
Q from most other popular models. We will calculate the high energy e + e~ spectrum in the quark-antiquark (q — q) 
center of momentum (CM) frame. Here, q stands for u, d, s, etc. Then, the result is convoluted to the q and q structure 
functions in the proton. 

For the q — q annihilation, the CM frame kinematics of the q — q pair is useful, 

q + q — > E-(E,p) + E+(E,-p) (28) 
— ► {e-(ei,qi) + Ni(E - ei,p - qi)} + {e + (e 2 , q 2 ) + N 2 (E - e 2 , -p - q 2 )} . 
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The q — q partem contribution in the proton-proton collision is given by 

a( PP ^e-e + X) = ^^J d ^J dx 2 f q {x x , Q)f s (x 2 , Q)^y^. (29) 

Here f q (xi,Q) and fq(x 2 ,Q) are parton distribution function(PDF) of q and q with p q = xip p and p g = x 2 p p , 
repectively. We use MRST99 C++ code for PDF[H]. Sum is for all quarks except for t and is neglected from now 
on. Q is relevent energy scale of PDF. s = x\x 2 s with y/s =14 TeV in the designed LHC energy. Aq is obtained by 
cross section at the parton level, and is given by 



A a = {Q q QeY + Q q QeRe(r)(sl + 4)(4 + e R ) + |rf[(4) 2 + (4)11(4)' + (4)1, 



(30) 



where Q q = +2/3 for up- type quarks Q q = —1/3 for down type quarks, Q e = — 1 are electromagnetic charge of up 
and electron, and 4 = 1 — | sin 2 6w, e q R = — | sin 2 6w for up-type quarks, e\ = — 1 + | sin 2 0^, 4 = +| sin 2 6*vk fo r 
down type quarks, 4 = — 1 + 2 sin 2 9w and e e R = 2 sin 2 #vf are Z boson couplings of quarks, and electron, respectively. 
6w is the weak mixing angle [HJ . r is given by 



V2G F Mf 



s - M% + iM z T z 47ra' 

where Gp is Fermi constant, Mz and are Z boson mass and decay rate. Eq. (|29p can be rewritten as 
a(pp — > e~e + X) = C [ ds [ dx\ [ dx 2 6(s — Xix 2 s)f q (xi, Q)f q (x 2 , Q)^r~ 



(31) 



fei / efo 2 — (5(.t 2 —)fq(xi,Q)fq(x2,Q)^J- 

X\S X\S s 

C I ds I d Xl — f g ( Xl ,Q)f q ( — ,Q)^M 



(32) 



I d ^~dl^ pp 



e~e + X) 



where 



f/lT 

rf.s 



(pp-f e~e+X) = C 



dx 1 —f q ( Xl ,Q)fq(—,Q)^-.. 



(33) 



which is cross section between s + ds and s. Here C = 87ra 2 /9. 

Similarly, ^(pp — > E~E + X) is obtained by simple calculations, which is given by 



1 



4,(5) 



daii f q (xi,Q)fq{ ,Q) 

s / s XlS XlS s 



!--2 

7 



3/2 



(34) 



where 7 = Epjmp. A' is given by 

4, = (Q q Q E ) 2 + Q,Q £ Jte(r)( e « + 4)(ef + e E R ) + |r| 2 [(4) 2 + (4) 2 ][( £ f ) 2 + {e E R )\ (35) 

where Qe = — 1 and 4 — e j? = 2 sin 2 6yy- The production cross section of E~E + is approximately a half of that of 
electron-positron pair at weak scale. Cross section is shown in Fig. [T] 

However, since E decays inside the beam line, electron-positron excess appears at lower energy scale. Therefore we 
can detect electron-positron pairs in collider. Cross sectrion can be rewritten in elctron-positron transverse invariant 
mass m ee = \J (qi + q 2 ) 2 = yeiT^T cosh(A?7) — qi^ • <\ 2 t where e^T, <liT and Arj are transverse energy, transverse 
momentum and rapidity difference, respectively. It is given by 



/d 
dm ee —^-(pp -► E-E+X 
dm ee 



where 



da 

dm et 



(pp -> E~E + X -> e~e + X) = a(pp -> E~E + X) 



dfl^dfl- 



dm e 



1 rfr_ 



-{E- 



iV) 



(36) 



(37) 
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FIG. 1: The EE production cross section. The horizontal axis is the CM energy in GeV. The unit of the vertical axis is 
GeV" 2 . 



Also, these events are accompanied by large missing energy and this physics is very similar to slcpton pair production 
of odinary MSSM. This is shown in Fig. fj 

The main difference of our calculation from the ordinary MSSM result is the produced £"s totaly decay to electrons 
but not to the other leptons. Due to this feature, we can use transverse momentum variable, tot2 to determine the 
mass of E (33[. Note that the phenomenology of E + E~ production is very similar to that of the selectron-antiselectron 
production because in both cases the heavy particle decays to e ± plus missing energy. So, it is crucial to distinguish 
them. But, the only way to distinguish them is to measure their spin at the LHC. Below, we concentrate on the total 
cross section but not the angular distribution. In fact, there are many analyses for measuring the particle spin at 
the LHC machine [13, [3f| . Since our E particle is the fermion, their polarization property looks the same as the KK 
fermion in the UED models, and it may be difficult to distinguish our fermion from the UED fcrmions. But we can 
expect to distinguish E from the selectron. For the same E and selectron mass, their kinematics are similar but angular 
distribution of electrons might be quite different in the two cases. If both E + E~ and the selectron-antiselectron pair 
are produced, we can observe two mj- 2 s which might be spectacular. 



V. CONCLUSION 

The recent PAMELA satellite data from the galactic source suggested the positron excess @ but no anti-proton 
excess [l3j . Regarding these observations, two dark matter components, the neutralino x an d the singlet fermion 
N, were suggested in Rcf. Q- The needed coupling has been given by Eq. (Q]). Here, we studied the flipped-SU(5) 
where charged SU(2) singlet leptons such as en, En and E R of Eq. (|TJ) can be its allowed representations, but a 
similar coupling for an anti-proton excess is not allowed because quarks must be embedded in SU(5) non-singlet 
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FIG. 2: The expected e + e~ production through E + E~ decay. The horizontal axis is the invariant pass of e + e~ pair in GeV. 
The unit of the vertical axis is GeV -2 . 



representations . 

In addition, there exists an interesting suggestion of the possible existence of VWLP with the coupling strength (J6j) 
from the study of white dwarf cooling [4[. In this paper, we tried to interpret this VWLP as an electrophilic axion with 
the axion decay constant F a ~ 1.4 x 10 10 GeV. In the DFSZ model, the corresponding F a would be ~ 1.2 x 10 9 GeV 
which would be barely consistent with the SN1987A bound. This kind of electrophilic axion can arise in 'variant very 
light axion' models where different families possess different PQ quantum numbers T. We argued that the inverted 
mass pattern of the first family from that of the second and third families is the logic for assigning different PQ 
quantum numbers for different family members. Allowing an effective domain wall number of ^ [See Footnote [TJ], 
we constructed an electrophilic axion. The symmetries we introduced in (J3J) is bigger than SU(5) xD(l)x xU(l)r, to 
allow a GUT scale VEV with an intermediate scale F a and R-parity. With these symmetries, we explained the mass 
hierarchy between the top and bottom quarks. In this model, the gauge couplings and the top quark Yukawa coupling 
are the only renormalizablc couplings. 

This model needs a charged SU(2) singlet lepton E at the electroweak scale. So, the existence of E is absolutely 
needed for this suggestion of the positron excess to work. Here, the discovery potential of E, at the LHC experiments 
by observing the electron and positron spectrum, is presented. We calculated the e + e~ spectrum at the center of 
momentum frame of the quark and anti-quark pair and compared it with the expected background. We pointed out 
that the SU(2) singlet E with the coupling JT]) is distinguishable from other possibilities we can imagine. Therefore, 
the discovery of E at the LHC experiments will indirectly confirm the two dark matter component hypothesis of • 
Or it can be ruled out from the LHC experiments by the absence of the e + e~ excess above the background. 
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